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Abstract 

In this paper, we introduce a generalized composite fading distribution (termed extended generalized- 
K (EGK)) to model the envelope and the power of the received signal in millimeter wave (60 GHz or 
above) and free-space optical channels. We obtain the first and the second-order statistics of the received 
signal envelope characterized by the EGK composite fading distribution. In particular, expressions 
for probability density function, cumulative distribution function, level crossing rate and average fade 
duration, and fractional moments are derived. In addition performance measures such as amount of 
fading, average bit error probabiUty, outage probabiUty, average capacity, and outage capacity are 
offered in closed-form. Selected numerical and computer simulation examples validate the accuracy 
of the presented mathematical analysis. 

Index Terms 

Composite fading distribution, generalized-K distribution, probability density function, cumulative 
distribution function, fractional moments, level crossing rate, amount of fade duration, moments, amount 
of fading, average bit error probability, average capacity. 

I. Introduction 

Radio wave propagation in wireless millimeter wave (60 GHz or above) and free-space optical 
(FSO) channels is a complicated phenomenon characterized by three interrelated phenomena: 
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path-loss with distance, shadowing (or long-term fading), and multipath (or short-term) fading. In 
virtue of fact that shadowing and multipath fading depend on reflection, refraction and scattering 
while the path-loss is distance dependent [1], [2, references therein], Gaussian-based models are 
usually employed to model these kind of channels due to their mathematical tractability [1], [2]. 
Indeed, the Gaussian approximation for the in-phase and quadrature components of the received 
random vectors leads to the commonly used channel fading models under short-term and long- 
term fading conditions, for the desired as well as the interfering signals. For instance, short- 
term fading distribution models include the well-known Rayleigh, WeibuU, Rice, Nakagami-m, 
generalized Nakagami-m models [3]-[6]. On the other hand, long term fading phenomena are 
modeled by the well-known lognormal distribution [7], [8]. These fading models are typically 
used to fit the histogram of the empirical/experimental measurements of the envelope of the 
received random signals. However, there exist fading situations, for which no distribution seems 
to adequately fit the experimental data, although one or another may yield a moderate fitting. This 
is accentuated for millimeter wave and FSO channels since multipath (small scale) fading and 
shadowing (large scale) fluctuations occur simultaneously leading to composite fading channels. 
In this sense, composite fading models are particularly important for the design of future wireless 
communications systems [9]. 

Modeling of composite fading channels is important to comprehend, contemplate and analyze 
several physical problems in wireless communications including interference effects in cellular 
systems, multiple input multiple output (MIMO) network, distributed antenna systems, coop- 
erative / multihop relay networks, and optical communications. In a typical signal propagation 
scenario, the received signal will show fading consisting of very rapid fluctuations around the 
mean signal level superimposed on relatively slow variations of the mean level. It is therefore 
not a misstep to assume that the local mean power of the multipath fading is a random variable 
(RV) distributed over (0, oo) according to a lognormal or gamma distribution. For example, in 
60 GHz non-line-of-sight (NLOS) propagation, the standard deviation of shadowing is typically 
larger than that of propagation at 5 GHz [10], [11] due to the fact that human body (moving 
objects) shadowing^ is a significant propagation effect in millimeter wave environment [12], [13]. 
However, the diffraction effects (propagation of the electromagnetic field behind the obstacle) 

'Reported in the literature that signal power can be considerable lost (i.e., up to 20 dB) if the line-of-sight (LOS) component 
is blocked by a person [12]. 
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are significantly smaller in comparison with 5 GHz band (shadowing zones are very sharp in 
60 GHz) [14]. In short, the environment changes very fast with respect to the diversity and 
average power of the received signal. In addition, in FSO communication systems, atmospheric 
turbulence, which originates from variations in the refractive index of the transmission medium 
due to inhomogeneities in temperature and pressure changes, increases the standard deviation of 
shadowing on the channel fading the received signal is subjected to. As a result, the identification 
of a tractable probability density function (PDF) to describe the physical problems at the 
background of shadowing effects is important in order to study the performance characteristics 
of emerging and future wireless communication systems operating in this kind of environments. 

In this paper, we focus on a new probability distribution, which is termed extended generalized- 
K (EGK), to model the fading in wireless millimeter wave channels and FSO environments. The 
EGK distribution has five parameters and is the extension of the composite fading distribution 
proposed in [15]-[18]. In addition, the EGK distribution has some good tail properties and 
includes most of the well-known fading distributions in the literature as either special or limit 
cases. 

The remainder of this paper is organized as follows. In Section II, the EGK distribution is 
introduced and some of its special cases, which are commonly used in the literature, are outlined. 
Moreover, the first order statistics of the EGK distribution are derived in closed-forms. Relying 
upon these first order statistics. Section in contains the second order statistics of the EGK 
distribution such as level crossing rate (LCR) and average fade duration (AFD). In Section IV, 
well-known performance measures for digital communications systems such as amount of fading 
(AoF), average bit error probabilities (ABEP), outage probability (OP), average capacity (AC) 
and outage capacity (OC) are derived for EGK fading channels, utilizing the results obtained in 
previous sections. Finally, conclusions are drawn in the last section. 

II. Extended Generalized-K Fading Channels 
Let the extended generalized-K (EGK) random process be defined as 

' -V ' 

= S{t) = X{t) 
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where Se{t) and X^{t) are two sets of independent Rayleigh distributed RVs with the average 
unit powers E [S'|(t)] = E = 1 such that S{t) and X{t) are a pair of independent 

RVs representing the shadowing (large-scale) and multipath fading (small-scale) components, 
respectively. Note that both S{t) and X{t) are distributed according to generalized Nakagami-m 
PDFs with the average powers E[S'^(t)] = ilg and E[X2(i)] = O^, where E[-] denotes the 
expectation operator. In (1), m (0.5 < m < oo) and ^ (0 < ^ < oo) represent the fading 
figure (diversity severity / order) and the fading shaping factor, respectively, while (0.5 < 
rrig < oo) and (0 < ^ < oo) represent the shadowing severity and the shadowing shaping 
factor (inhomogeneity), respectively. In addition, the parameters (3 and are defined as /3 = 
r (m + 1/^) /r (m) and = F (m^ -|- 1/^^) /T (m^), respectively, where T (•) is the Gamma 
function [19, Eq. (6.5.3)]. 

A. Probability Density Function 

The PDF of the received signal envelope TZ introduced in (1) is a product of two RVs (that 
is, TZ = SX) and as such, using the Mellin transform [20], [21], it is represented in a general 
compact form given in the following definition. 

Definition 1 (Extended Generalized-K RV). The distribution TZ follows an extended generalized- 
K (EGK) distribution if the PDF ofTZ is given by 

where the parameters m (0.5 < < oo), ^ {0 < d < oo)- ''^s (0-5 < < oo), (0 < 
^f, < oo) and Q = ^g^x (0 < < oo) are defined above, and where T (-, -, -, ■) is the 
extended incomplete Gamma function defined as T (a, x, 6, /3) = r"~^ exp (— r — hr~^^ dr, 
where a, /3,b e C and x e R+ [22, Eq. (6.2)]. 

In what follows, we utilize a shorthand notation, i.e., IZ ^ KLg (m, ^, m^, Vt) which denotes 
that TZ follows a EGK distribution with the fading figure m, the shaping factor ^, the shadowing 
figure nig, the shadowing shaping factor and the average powers Jl. 

To the best of our knowledge, the EGK distribution 7?. ~ /Cg (m, ^, m^, Jl) has the advantage 
of modeling the envelope statistics of most known wireless and optical communication channels. 
For example, by using [22, Eq. (6.41)] and setting ^ = 1 and ^ = 1 in 7?. ~ /Cg (m, nig, Vl), 
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we get the PDF of generalized K (GK) composite distribution, namely [16] 



Pn{r) 



4 /mrrig 



) 



m+m, 
2~' 




(3) 



r(m)r(mj V Q 



m—m, 



Moreover, as shadowing severity approaches infinity (m^ — > oo), which means that the 
distribution of the shadowing S{t) part in (1) follows the dirac's delta distribution (ps{r) — 
S {r — Q)), the PDF (2) simplifies into the PDF of the generalized Nakagami-m^, given by [6] 



whose the special or limiting cases are well-known in literature as the Rayleigh (m^ = 1, = 1), 
exponential (m^ — 1,^^ — 1/2), Half-Normal (m^ = 1/2, = 1), Nakagami-m (^^ = 1), Gamma 
(^^ = 1/2), WeibuU (m^ ==1). lognormal (m^ ^ oo,^^ ^ 0), and AWGN (m^ oo,^^ = 1). 
For the other commonly used channel fading models, special or limiting cases of the EGK 
distribution are listed in Table 1. Regarding this great versatility of the EGK distribution, it is 
important to notice that the EGK distribution offers a kind of unified theory to statistically model 
the envelope statistics of most known wireless/optical communication channels. 

Note that all distributions listed in Table I as the special or limiting cases of the EGK 
distribution are proposed purely from empirical fitting of measured data to a statistical dis- 
tribution with their corresponding tail properties lim P7^ (r) and lim dp-ji{r)/dr. If the tail 
properties approximate to zero for such high r values, i.e., lim,.^oo Pt^ = 0+ > and 
lim dpTz{r)/dr = 0~ < 0, then the probability of low amplitude values increase since the 



condition pTi{r)dr = 1 is always valid. Explicitly, these tail properties emphasize the 
distribution at low amplitude values as seen in Fig. 1. For example, the variation of the EGK 
PDF pti (r) is depicted in Fig. 1 for 3 dB thresholds (that is, r = V2Q, r — and r = 
in order to accentuate how the tails of EGK PDF p-n (r) curves changes with respect not only 
to both fading figure m and fading shaping factor ^ but also to both fading shaping factor 
^ and shadowing shaping factor As seen in Fig. 1(a), as either fading figure m or fading 
shaping factor ^ goes to low possible values (m ^ | or ^ ^ 0), the tails of EGK PDF 
Pn (f) goes to zero, which means that the probability of the distribution at low amplitude values 
increases since Pn {r)dr = 1 as mentioned before, and the performance measures of wireless 

^Note that the generalized Nakagami-m RV is the square root of the generalized gamma RV proposed by Stacy [6]. 
December 12, 2010 DRAFT 
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communications systems consequently deteriorates. Additionally, the same consequences are 
analogously extrapolated for both fading shaping factor ^ and shadowing factor as seen in 
Fig. Kb). 

In additive white Gaussian noise (AWGN) channel, the distribution of the instantaneous signal- 
to-noise ratio (SNR), Q — TZ^/Nq can be directly expressed in terms of average SNR, 7 = 
E[^?] = E[7^2]/A^o = Q/Nq and A^o representing the power of AWGN noise. The PDF of 
Q, which we term the extended generalized gamma (EGG) PDF, is defined in the following 
definition. 

Definition 2 (Extended Generalized Gamma RV). The distribution Q follows an extended gene- 
ralized gamma (EGG) distribution if the PDF of Q is given by 

Pg{l) = z- — TT^r^ — 7^ r -m-,0, — 7% ^ • (5) 

r(mjr(m) V7y \ is \ 1 J is J 

In what follows, we utilize a shorthand notation, i.e., Q ~ {rn,^,mg,^g,^) which denotes 
that 7 follows a EGG distribution with the fading figure m, the shaping factor the shadowing 
figure nig, the shadowing shaping factor and the average powers 7. 

B. Fractional Moments 

The fractional moments E [TV'] , k e 1R+ are crucial for several reasons. First, purely moment- 
based measures, such as the average SNR and AoF [23], which can be computed using only the 
first and second central moments of the SNR at the diversity combiner output, are commonly used 
to characterize the diversity systems. In addition, as shown in [24], higher order moments can 
also be used to characterize the statistical behavior of the output instantaneous SNR distribution 
for certain diversity systems. Secondly, more widely used performance measures such as ABEP 
(which is suitable for digital modulations) and OP, can be typically computed as shown in 
[25] using the Laguerre moments (computed based on fractional moments) in the case of that 
the moment generating function (MGF) and PDF of the received instantaneous SNR are not 
available. 

The fractional moments E [TZ^] , k e M"*" of the signal envelope TZ ^ ICg (m, m^, fl) are 
given in the following theorem. 
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Theorem 1 (Fractional Moments of EGK RV). For k e 1R+, the kth moment E [U''] of the 
signal envelope TZ ~ ICg (m, ^, m^, C,s, ^) is given by 



(6) 



rK)r(m) 

Proof: Utilizing [22, Eq. (6.22)] with [26, Eqs. (2.1.4), (2.1.5) and (2.1.11)], the PDF ptz (r) 
can also be represented as 

2 



Pn[r) 



r(mjr(m)r 



H 



2,0 
0,2 



(7) 



where H^^"[-] is the Fox's H function^''*, and where means that the parameters are absent. 

Using (7), the kth moment E — r'^pn (r) dr can be written as 



' ^0,2 



dr. 



(8) 



From the definition of Mellin transform [26, Eq. (2.5.1)] and using [26, Theorem 2.2], the /cth 
moment E \JZ^'\ of the signal envelope TZ Kg (m, ^, m^, Jl) is readily obtained as in (6), 
which proves Theorem 1. ■ 

Let us consider some special cases of (6) in order to check its analytical correctness. When 
setting the fading shaping factor ^ = 1 and the shadowing shaping factor = 1, (6) reduces into 
[16, Eq. (7)] as expected. For the other commonly used channel fading models, which are listed 
in Table I as the special or limiting cases of the 7?. ~ /Cg (m, ^, m^, O), the k\h moment can be 
readily obtained; for example, when the shadowing severity approaches infinity (m^ — > oo), 
(6) simplifies into the k\h moment of the generalized Nakagami-m RV [29, Eq. (5)]. 

Using Theorem 1, the k\h fractional moment E \Q'^^ , k e Mf^ of Q Qg (m, ^, m^, 7) can 
be readily obtained as 



E [g^] = E [n^''] 



r(m, + A)r(m+|) 
rK)r(m) 



7 



(9) 



^For more information about the Fox's H function, the readers are referred to [26], [27] 

"'Using [28, Eq. (8.3.22)], the Fox's H function can be represented in terms of the Meijer's G function [28, Eq. (8.2.1)] which 
is a built-in function in the most popular mathematical software packages such as MATHEMATICA®. 
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C. Cumulative Distribution Function 

The CDF P-ji (r) of the signal envelope TZ is defined as the probability that the received 
envelope TZ falls below a threshold level r of the signal envelope. In the following theorem, the 
CDF ofRr^lCg (m, ^, m^, Vi) is given. 

Theorem 2 (Cumulative Distribution Function of EGK RV). The CDF P-ji (r) of the signal 
envelope TZ ~ Kg {171,^,111^,^^,^) is given by 



Pn{r) 



rK)r(m) 



tt3,1 
^1,3 



(1,1) 

K,^),(m,|),(0,l) 



(10) 



Proof: Upon having (7), i.e., the PDF of TZ ^ ICg (m, ^, m^, O) in terms of Fox's H 
function, and utilizing the equality [20, Eq. (4.18)], one can readily obtain the CDF P-n (r) = 
JqPv, (r) du of TZ JCg (m, ^, m^, Q) as in (10), which proves Theorem 2. ■ 
Note that the PDF of 7?. ~ /Cg (m, ^, m^, Q) given by (10) can also be represented as 



Pnir)^l- 



r(mjr(m; 



tt3,0 
"lll,3 



(1,1) 

(m„^),(m,i),(0,l) 



(11) 



by means of using the equality given in [20, Eq. (4.17)]. ExpUcitly, both (10) and (11) are 
identical functions, but represented differently in terms of Fox's H function. Let us consider 

special cases of (10); for example, when the shadowing severity approaches infinity (m^ 
00), we get [29, Eq. (3)] as expected, by mean of substituting lim p^"^''K ^ 1, where |6| <C a 
and |c| ^ a, into the Mellin-Bames integral representation [26, Eq. (1.1.1)] of the CDF given 
by (10) (or (11)). 

Note that, in the case of the difficulty in the computation of Fox'H function, (10) can embody 
the Meijer's G representation for the rational values of the parameters ^ and (that is, we 

let ^ = k/i and = kg/t^, where k, kg and are arbitrary positive integers.) through the 
medium of algebraic manipulations utilizing [28, Eq. (8.3.2.22)], namely 



Pn{r) 



r(mjr(m) 



G 



2\ksk 



^(-ksk) 



(12) 



with $ = y(2^i^2:^^^^+w (;t4)rn.-i(;(.^£)m-i and * = (M)^(^4)fe, where G^f[-] is the Mei- 



jer's G function [28, Eq. (8.3.22)], and the coefficients are defined as = ^, ^±1, 



a+m—l 



with a e C and m e N. Additionally, one can also readily approximate the CDF P-ji (r) with high 
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accuracy as the sum of PDFs p-ji (r) by means of employing GCQ formula [19, Eq. (25.4.39)], 
which converges rapidly and steadily, requiring few terms for an accurate result. Then, the definite 
integral in definition of CDF P-ji (r) = Pn (r) du can be accurately estimated as 



Pn{r) 



2m^ 



n=l 



(13) 



where ^ sin (^f^7r) and 0„ = | + | cos (^|^7r) for n e {1,2,..., N}, and where the 

summation stop index N could be chosen as A?^ = 30 or more to obtain a high level of accuracy. 

Again, using Theorem 2, the CDF Pg (7) of Q Qg (m, ^, m^, 7) can be readily obtained 
by substituting (10) into Pg (7) = Pn {^) and then changing ^ 7, that is. 



rK)r(m) 



tt3,1 



7 



-7 



(1,1) 

(m„^),(m,i),(0,l) 



(14) 



D. Moment Generating Function 

Upon utilizing both alternative exponential forms (i.e., Craig formula forms) of the Gaussian 

error function, i.e., erfc {^/x) = (2/7r) JJ"^^ exp (— xcsc^ {6))d6 [5, Eq. (4.2)] and crfc^ (v^) = 
(4/7r) /J'^'' exp (-xcsc^ (^))(i^ [5, Eq. (4.9)], the moment generating function (MGF) of ^ ~ 
Qg (m, ^, rrig, 7) is required regarding to the average symbol error probabilities of the receivers 
operating in wireless communications channels. The MGF of Q ^ Qg {m,^,mg,^g,^) is given 
in closed-form in the following theorem. 

Theorem 3 (Moment Generating Function of EGG RV). The CDF M.g{s) of the instantaneous 
SNR Q Qg (m, ^, m^, 4, 7) is given by 



Mg {s) 



1 



r(m^)r(m) 



tt2,1 
^^1,2 



7s 



(1,1) 
(m3,^),(m,l) 



(15) 



Proof: Upon having (7), i.e., the PDF of ^ ~ ("^, ^s, C 7) terms of Fox's H 
function, and utilizing the equality [20, Eq. (3.8)], one can readily obtain the MGF Mg (r) = 
exp (—57) Pg (7) (^7 of ^ ~ Qg (m, ^, m^, 7) as in (15), which proves Theorem 3. ■ 
Note that, upon following the same steps in the derivation of (12), (15) can be represented 
in terms of Meijer's G function, by favor of [28, Eq. (8.3.2.22)], for the rational values of the 
shaping parameters ^ = k/i and ^^/c^/^^, where k, £, kg and £s are arbitrary positive integers. 
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that is. 



G 



r(mjr(m) 



7* 



2\ksk 



^{-k,k) 

r^ilsk) .^{kst) 



(16) 



In addition, setting fading and shadowing shaping parameters to one (i.e., = 1 and = 1) in 
both (15) and (16), and then using [28, Eq.(8. 3.2.21)] produces the MGF of the instantaneous 
SNR at the output of a receiver operating over generalized-K channels, namely 

1 



Mg (s) = 



p2,l 



rrigm 


1 


7S 





(17) 



r(mjr(m) 

Regarding the special or limiting cases of the EGK distribution listed in Table 1, the MGF 
of the other commonly used channel fading models can be readily found by substituting the 
corresponding values of the fading figure m, the shaping factor ^, the shadowing figure m^, the 
shadowing shaping factor into either (15) or (16). 

III. Second Order Statistics 

Referring to (1), when S{t) and X{t) are signal envelopes in some scattering fading channel 
subjected to Doppler effect in wireless fading channels or to the turbulence effect originate from 
variations in the refractive index of the transmission channel in free-space optical channels, the 
signal envelopes S{t) and X{t) are time-correlated random processes and their time derivatives 
S{t) and X{t) are not independent from their corresponding envelopes but following zero-mean 
conditional Gaussian PDFs with respective conditional variances 



^2 



'x\x 2e V /5 



(18a) 



(18b) 



where a^^ = u:vsj\pl is with the relative vehicle speed vs with respect to the obstructs / keyholes 
causing shadowing and the wave number given by cu = 27r/A in terms of the wavelength A such 
that ujvs = 271 fs where fu is the maximum Doppler frequency shift influencing the shadowing 
component. Moreover, in (18b), cr^^ = ujux/V2 is with the relative vehicle speed px with 
respect to the transmitter such that uux — Stt/x where fx is the maximum Doppler frequency 
shift influencing the multipath fading component. In the following, identical mean power shall 



be assumed at all components (i.e., a 



•S'fe 



al, and a% 



a% for k ^ £). The time derivative of 

A.- 
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the received envelope TZ Kg (m, ^, m^, Q), i.e., TZ can be obtained by TZ — SX + SX — 
STZ/S + SX. Note that the time derivative 7t is a zero mean Gaussian RV with the conditional 

%al^^ + S^al.^^, where ct|,^^ = ct|,^ and cr^ = cr^. ,^ = a^., ^. For 



variance a. 



2 



fixed S = u and 71 = r, the conditional variance is readily obtained as 

9 \«= \« . .2e(l/e-l) 

Remembering that the time derivative 7^ is a zero mean Gaussian RV with the conditional 
variance given in (19), the conditional PDF ^7^17^5 {f\7lS) can be readily given by 

1 ( P \ 

Pn\ns {r\nS) = -= exp -— (20) 

distributed over — cxd < r < cxo. 

Finally, utilizing the conditional PDF (20), the second order statistics (i.e., the LCR and the 
AFD) of the fading envelope 7?. ~ /Cg (m, ^, m^, Q) can be obtained as we show in what 
follows. 

A. Level Crossing Rate 

The LCR of the fading envelope IZ ^ K,g (m, ^, m^, fi) at threshold r is defined as the 
rate at which the EGK random process crosses level r in the negative direction. Using the joint 
PDF of the EGK envelope random process TZ and its time derivative iZ, i.e., p-j^-f^ (r, r), we can 
obtain the LCR as follows 

poo 

Ln (r) = / rp^^n ('^' ^) ^r, (21) 







where the joint PDF p^^^^ (r, r) can be expressed in terms of conditional PDFs as follows 

POD 

Pnn ^) = / PnVRS Pn\s{r\u)ps{u)du, (22) 

Jo 

where ^7^17^5 (f|r, m) is the conditional Gaussian PDF as mentioned before. Moreover, ps{u) is 
the generalized Nakagami-m PDF representing the distribution of the shadowing component of 
the received envelope while pnisifW) = px ('^/^) /u. Substituting (22) into (21) and using 
After some algebraic manipulations, we have 

1 /"°° 

Liz (r) = -j= aT^^jisPTzis{r\u)ps{u)du. (23) 
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Substituting (19) into (23), we obtain the exact solution for the LCR as in 




Let us consider some special cases in order to check analytical simplicity and accuracy of 
our result. As seen in Table I, setting the shadowing and multipath fading shaping factors 
(inhomogeneity) ^ = 1 and = 1, respectively, in (24) results in the LCR of generalized-K [30, 
Eq. (9)], as expected. It is useful to mention that the above integral can be computed numerically 
with desired accuracy (e.g. by using the most popular mathematical software packages such as 
MATLAB®, MATHEMATICA® and MAPLEtm.), or alternatively, it can be readily estimated 
accurately in virtue of the GCQ rule [19, Eq. (25.4.39)] after changing the variable of the integra- 
tion in (24) as M — > tan(u). In addition, the above integral can be approximated by means of the 
Laplace approximation [30], [31]. Distinctively, we obtained the above integral by utilizing the 
Taylor series expansion y/l + x = limiv^oo Xl^Lo Cn) (^'"^"^ {1 — x) + x^~"'9 (x — where 
6 (•) is the Heaviside's theta (unit) function [32, Eq.(1.8.3)]. Shortly, applying the Taylor series 
expansion of + x and performing some algebraic manipulations, (24) can be converted into 
the sum of two incomplete integrals converging very fast, that is. 




where the auxiliary function Q„ (r) is defined as 
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Eventually, substituting [22, Eqs.(6.1) and (6.2)] into (25) results the closed form of the LCR 
for the fading envelope 7?. ~ /Cg (m, ^, m^, il) as follows 

e-„ (r)7 fp-„ - e M, 3^ Qi w, i) + 

W r f T'n - (r), (r), pi ] (27) 

with Vn = {jn^ — n) — ■f{m + n), where 7 (■,■,■,■) is the extended (lower) incomplete Gamma 
function defined as 7 (a, x, h, (3) = r"~^ exp (— r — br^^) dr, where a, /3, 6 G C and a; G M+ 
[22, Eq. (6.1)]. 

In Fig. 2, the LCR for EGK fading channels is depicted for different parameters. The shad- 
owing increases when the shadowing figure ms — > | and the shadowing factor oo> or it 
decreases when the shadowing figure ms — > 00 and the shadowing factor — > 0. As seen in 
Fig. 2, the worst case shadowing is represented by ms — > | and = 1. The shadowing and 
multipath fading figures {771^,771} and/or {^^,0 increase, the LCR decreases as expected. 

It may be useful to note that the LCR given in (27) converges rapidly. Indeed, choice N = 1 
is enough for an approximation. As such, (27) can be readily approximated as 

\j^0Ai r(m)r(m.) \ 

S-,Wr(n-«ii.S^M.||Q(.i)i) + 

ic.(r)r(p,-«il.S^(,.).||c.wi)}. (28) 
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B. Average Fade Duration 

The AFD of the received envelope IZ ~ ICg (m, ^, m^, Vt) at threshold r is defined as the 
average duration that the EGK random process representing the envelope of the received signal 
remains below the level r after crossing that level in the downward direction, that is, 

T{r) = ^^y (29) 

Finally, substituting (10) and (27) into (29), one can readily obtain the AFD of the received 
envelope TZ ~ /Cg (m, m^, f2) at threshold r. The AFD of the received envelope TZ ~ 
Kg (m, ^, rrig, Vt) at threshold r is depicted in Fig. 3 for some scenarios. 

IV. Performance Measures over EGK Fading Channels 

In this section, using the expressions obtained in previous sections, the expressions of AoF, 
ABEP, OP, AC, and OC for digital communications systems operating over EGK fading channels 
are given in closed-forms, and their analytical simplicity and accuracy are checked by simulations. 

A. Amount of Fading 

The AoF introduced in [23] as a unified performance measure of the severity of fading, is 
an important measure for the performance of a diversity system since it can be utilized to 

parameterize the distribution of the signal-to-noise ratio (SNR) of the received signal. Referring 
to (9), the AoF of Q ^ Qg (m, ^, m^, 7) can be computed by 

var [g] _ E|^2i 

Then, the AoF of Q ^ Qg (m, m^, 7) is given in the following theorem 



AoF = = - 1 (30) 



Theorem 4 (Amount of Fading of EGG RV). The AoF of Q Qg {m,^,mg,^g,^) is given by 

r (mj r fm, + r (m) T (m + j) 
AoF = } ^-^ ^ ^ - 1. (31) 

Proof: Using (6) , the proof is obvious. ■ 
As readily seen in (31), when the fading figures {mg,m} and/or {C^,^} increase, the AoF 
decreases as expected. 
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B. Average Bit Error Probabilities 

The instantaneous bit error probabilities (IBEP), conditioned on the instantaneous SNR Q ~ 
Gq T^si Cs) 7)> ill ^ AWGN channel may be written in compact form as [5, Eq. (8.100)] 



(32) 



where a depends on type of modulation scheme (| for orthogonal FSK, 1 for antipodal PSK), 
h depends on type of detection technique (^ for coherent, 1 for non-coherent), and V (•, ■)is 
moreover complementary incomplete Gamma function [19, Eq. (6.5.3)], respectively. Using the 
PDF of the instantaneous SNR ^ ~ (m, ^, m^, 4, 7), i.e., (5) and utilizing Pe {Q), the ABEP 
in EGK fading channels is given by 



'0 



(33) 



Eventually, substituting (32) and (5) into (34), the ABEP in EGK fading channels is given in 
closed-form in the following theorem. 

Theorem 5 (Average Bit Error Probabilities in EGK Fading Channels). Let Q ^ Qg (m, 7) 
be the instantaneous SNR at the output the receiver operating in wireless EGK fading channels. 
Then, the ABEP Pe of the receiver is given by 



1 



-H: 



2,2 



for wireless EGK fading channels. 



07 



(1-6,1) (1,1) 

(m„i),(m,i),(0,l) 



(34) 



Proof: Substituting (32) and (5) into (34), the ABEP Pe can be given as 



E 



7 



, 7"^^-'r(6,ar)r m,-m^,0, 
2r(6)r(mjr(m) Jo \ L \ 1 



7«,^M7 (35) 



where substituting Meijer's G and Fox's H function representations of both the incomplete gamma 
r (•, •) [28, Eq. (8.4.16.2)] and the extended incomplete gamma F (•, •, •, •) [22, Eq. (6.22)], we 
get 



1 



2r(6)r(mjr(m)7o 7 



n2.o 
"^1,2 



07 



1 

0,6 



H, 



2,0 
0,2 



7 



-7 



(m,, ^),(m,)) 



d-f (36) 
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Finally, representing the Meijer's G function in (36) in terms of Fox's H function representation 
[26, Eq. (2.9.1)] and applying [26, Theorem 2.9] onto (36) results as in (34), which proves 
Theorem 5. ■ 
To check the analytical correctness, let us consider some special cases of (34). Upon setting 
the fading shaping factor ^ ^ 1 and the shadowing shaping factor ^ 1 and utilizing [26, Eq. 
(2.9.1)], (34) is simplified to 

^ 1 



E 



-G 



2,2 



171^171 


1 - 6, 1 " 




rrig, m, 



(37) 



2r (6) r (mjr (m) ^-^ 

which is the ABEP Pe of the receiver operating in well-known wireless composite fading 
channel, i.e., wireless generalized-K fading channels. 

When the shadowing severity rrig approaches infinity (nig — > oo) meaning that there is no 
shadowing effect in the channels, (34) is simplified to the well-known ABEP of binary digital 
modulation schemes in generalized Gamma fading channels [33, Eq. (10)], namely 



E 



^ tt1,2 

2r(6)r(m) 



(1-M),(1,1) 
(m,|),(0,l) 



(38) 



Again, in the special case when the fading shaping factor ^ — > 1 and the shadowing severity 
rrig ^ oo , (34) becomes the more familiar expression for the ABEP in a flat fading Nakagami-m 
channel (see [5, Eq. (8.106)], for example), i.e., using [26, Eq. (2.9.1)], [26, Eq. (2.9.16)] and 
[28, Eq. (6.7)], we get 

1 



Pf 



2r (b) r (m) 



pl,2 





1-6,1" 


07 


m, 



r (/// + b) 



m 



1, m + 6; m + 1; 



m 



m + a'j 



(39a) 
, (39b) 



where 2^1 



m + J 2r (6) r (m + 1) \m + 07^ 

■] is the Gaussian hypergeometric function defined in [28, Eq. (7.2.1.1)]. Fur- 
^ = 1, (39a) further simplifies to the well-known result 



thermore, for the special case m 

b 



E 



[5]. 



2 \ m+a-y 

The analytical simplicity and accuracy of (34) are checked by simulations for the noncoherent 
BFSK (NCFSK)/ differentially coherent BPSK (DPSK) and coherent BPSK/BFSK in Fig. 4. As 
seen in Fig. 4, when the fading figure increases, then the ABEP decreases since increasing the 
fading figure decreases the fading severity of the channels. Furthermore, as the shape parameter 
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increases the PDF of the fading becomes skewed around the average power, which impUes that 
the ABEP decreases. 

C. Outage Probability & Outage Capacity 

The OP is defined as the probabiUty that the instantaneous error rate exceeds a specified value 
or equivalently that the instantaneous SNR Q falls below a certain specified threshold ^th, i-e.. 



where pg (7) is the PDF of the instantaneous SNR Q (see (5)). In other words, the OP Pout is 
the CDF of the RV G evaluated at ^th- Then, using (14), the OP Pout is easily obtained as Pout — 
Pg {'fth)- In Fig- 5, by using Pout — Pg {ith) as a function for the normalized outage threshold 
1th/ 1, Pout is depicted using different fading figure m, fading shaping factor ^, shadowing 
severity and shadowing shaping factor Clearly, as seen in Fig. 5, Pout decreases as m or 
^ increases because increasing m decreases the fading severity and increasing ^ skews the PDF 
of the fading around the average power 7. 

The OC is another important statistical measure to quantify the spectral efficiency in fading 
channels, which is defined as the probability that the instantaneous capacity, Cg falls below a 
certain specified threshold Cth, i.e., Cout = Pr (0 < < Cth)- With the aid of [34, Eq. (5.1)], 
the OC C^t can be given in terms of (14) as follows C^t = Pg (2'^*''/^ - l). 

D. Average Capacity 

We consider a signal's transmission of bandwidth W over AWGN channel, then the Shannon 
capacity is defined as Cg {Q) = VFlog2 (1 + Q), where Q ^ Qg (m, ^, m^, 7) is the received 
instantaneous SNR. Then, the AC Cg can be obtained by averaging Cg{Q) with the PDF of 

g^gg (m, 2e, m„ 2C, 7) with 7 = E [TZ^/No] = n/No, i.e.. 



^0 

where pg (7) is given in (5). With the aid of [22, Eq. (6.22)], utilizing the Fox's H function 
representation of (5) and substituting [20, Eq. (2.4.3/1)] into (41), we can obtain the AC in EGK 



Pout = Pr{0<g <^tt,) = / pg{l)d^ 



(40) 



Jo 




(41) 
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fading channels in the closed form as 

W 



Cc 



log (2)r(m)r(mJ 



H 



4,1 
2,4 



7 



(o,i: 

K,^),(m,i),(0,l),(0,l) 



(42) 



Note that, using [28, Eq. (8.3.21)], for the shaping factors C = C = 1' (42) simplifies to the 
average capacity of the generalized-K fading channels, namely 



W/\og{2) _^4,i 



r (m^) r (m) 



2,4 





0,1 


. 7 


nig, m, 0, 



(43) 



Note that the AC of the generalized-K fading channels given in (43) is in a more compact form 
than in [18, Eq. (9)]. In addition, when we set the shadowing shaping and figure parameters 
as — 1 rrig ^ oo, referring to Table I and using [19, Eq. (6.1.46)] in the Mellin-Bames 
representation of (42) [26, Eq. (1.1.1)], (42) simplifies into the AC of Nakagami-m fading channels 



given by C. 



G 



3,1 



T log(2)r(m) ^2,3 



rn 


0,1 




m,0,0 



[35, Eq. (3)]. As seen in Fig. 6, improves with an 



increase of fading figures {nig^m} and shaping factors {^5,^}, as expected. 

V. Conclusion 

In this paper, we introduced a very general composite fading distribution to model the envelope 
and power of the received signal in fading channels, which we term the EGK composite fading 
distribution. We also studied second order of statistics such as level crossing rate and average 
fade duration for EGK fading channels. We finally obtained closed-form expressions for the 
amount of fading, average bit error probability, outage probability, outage capacity and average 
capacity of the EGK fading channels. 
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(a) With respect to the fading severity m and the fading (b) With respect to the fading shaping factor ^ and the 
shaping factor 5. shadowing shaping factor 



Fig. 1. The variation of the EGK PDF p-n (r) for TZ ~ Kg (m, 5, m^, ^s, Q.) is given with respect not only to both fading figure 
m and fading shaping factor ^ but also to both fading shaping factor ^ and shadowing shaping factor at 3dB envelopes: (i) 
r = yin, (ii) r = (iii) r = ^/flJ2, where = E [7^^] . 
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Fig. 2. Normalized LCR versus normalized envelope level in EGK fading channels, where the doppler frequency shifts of 
multipath and shadowing components of the fading are assumed fd = fs = fx- 
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_14 -12 -10 -8 -6 -4 -2 2 4 



r/v^E]^ [dB] 

Fig. 3. Normalized AFD versus tire normalized envelope level r/-y/E [R,^] in EGK fading channels, where r is the level and 
E [TZ^] is such a average power as E [IZ^] = i^x^s- 
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Average Power (SNR), 7 



Fig. 4. Average bit error propbabilities of BPSK, DPSK and NCFSK binary modulation schemes over EGK fading channels, 
i.e., analysis of Eq.(34). 
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Average Power (SNR), 7 



Fig. 5. Outage probability over EGK fading channels, i.e., analysis of Eq.(40) using (14). 
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TABLE I 

Some Special Cases of the EGK Distribution. 
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